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(54) Iterative channel estimation in multicarrier receivers 



(57) A method of channel estimation comprises the 
steps of receiving a signal and applying a fast fourier 
transform (FFT) toobtain apluralityof frequency domain 
sub-carriers. Probabilities of coded bits are produced 
and channel coefficient estimation is performed for each 



of the plurality of frequency domain sub-carriers using 
channel coefficient estimates of other ones (preferably 
all others) of the plurality of frequency domain sub-car- 
riers, which are based on the probabilities on coded bits. 
An OFDM communication unit and OFDM communica- 
tion system are also provided. 
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Description 

Field of the Invention 

5 [0001 ] This invention relates to an estimation of radio-channel propagation conditions, in order to equalise their effect 
on transmitted data. The invention is applicable to improving an accuracy of channel estimation in an orthogonal fre- 
quency division multiplex (OFDM) wireless receiver. 

Background of the Invention 

w 

[0002] Wireless communication systems, for example cellular telephony or private mobile radio communication sys- 
tems, typically provide for radio telecommunication links to be arranged between a plurality of base transceiver stations 
(BTSs) and a plurality of subscriber units, often termed mobile stations (MSs). 

[0003] Wireless communication systems are distinguished over fixed communication systems, such as the public 
15 switched telephone network (PSTN), principally in that mobile stations move between BTS (and/or different service 
providers) and in doing so encounter varying radio propagation environments. 

[0004] Methods of communicating information simultaneously exist where communication resources in a communi- 
cation network are shared by a number of users. Such methods are termed multiple access techniques. A number of 
multiple access techniques exist, such as frequency division multiple access (FDMA), time division multiple access 
20 (TDMA), and code division multiple access (CDMA). 

[0005] The communication link from a BTS to a MS is referred to as the down-link. Conversely, the communication 
link from a MS to the BTS is referred to as the up-link. 

[0006] In wireless communication systems, there is a need to estimate the effect caused by the 'wireless communi- 
cation channel' on the data being transmitted. Channel estimation is required so that the received data can be equalised 
25 to reduce, restore or minimise signal degradation caused by such transmission channel impairments. A radio channel 
is said to be equalised if the channel impairments can be eliminated or significantly reduced. 

[0007] In the field of this invention, the most conventional strategy in estimating radio channel propagation conditions 
is to model the radio channel by a finite impulse response (FIR) filter. The channel estimation is typically performed by 
periodically transmitting a known data sequence, generally referred to as a training sequence, over the desired radio 
30 channel. Such training sequences are a-priori known by the receiver. The training sequence is extracted from the 
received, desired data stream and is used to compute channel estimates. 

[0008] In the field of this invention, namely that of Orthogonal Frequency Division Multiplexed (OFDM) systems, also 
termed multi-carrier systems, the data stream is divided into several (N) sub-streams. These sub-streams are trans- 
mitted on 'N' orthogonal sub-carriers at different frequencies, for example by means of an inverse fast Fourier transform 

35 (IFFT) at the transmitting unit. 

[0009] One of the main advantages of OFDM systems is their very simple equalization scheme, which is reduced to 
a multiplication of the FFT outputs (of each carrier frequency) by the frequency domain channel coefficients. Of course, 
the true channel coefficients are unknown, and have to be estimated. As an illustration, the mean square error between 
the true channel coefficients and the channel coefficients at time 0 (time 0 being chosen when a training symbol is 

40 received) is plotted in function of time on FIG. 1 , for an HIPERLAN/2 channel with a mobile terminal moving at 3m/s. 
[0010] Referring now to FIG. 2, a classical Bit interleaved and convolutionally coded (BICC) OFDM system 200 is 
shown. The BICC OFDM system 200 includes a BICC OFDM transmitter 21 0. The BICC transmitter 21 0 receives a 
data stream 212 d, and convolutionally codes 215 the data stream to produce an output b, 217. The convolutionally 
coded output bj 21 7 is input to a bit interleaver 220. The output from the bit interleaver 220 is then input to a mapping 

45 function 225, which associates a subset of bits x k 228 to a location on a constellation. The type of constellation generated 
by the transmitter 21 0 is dependent upon the modulation scheme employed by the OFDM modulator 230. The OFDM 
modulator then outputs the OFDM modulated signal over the communication channel 235. 

[0011] A BICC OFDM receiver 250 receives the OFDM modulated signal over the communication channel 235. In 
effect, the receiver 250 performs the inverse operations of the BICC OFDM transmitter 21 0. In this regard, the receiver 

50 250 includes an OFDM demodulator 255 to translate the received constellation locations into a sequence of bit subsets. 
The bit subsets are then input to a demapping function 260. In the demapping function, bit metrics are computed to 
feed the Viterbi decoder 270. This computation involves the frequency channel coefficients, which must be estimated 
periodically. A bit deinterleaving function 265 arranges the received bit-stream for decoding in the Viterbi decoder 270. 
[0012] The above receiver arrangement presents a number of significant drawbacks. In particular, the process to 

55 demap, decode and perform channel estimation separately, as three distinct and independent functions, only enables 
the respective individual operations to be optimised. 

[0013] In J. Boutros, C. Lamy, F. Boixadera, "Bit-interleaved coded modulation for multiple-input multiple-output 
channels", IEEE ISSSTA , Sept. 2000, a method performing joint de-mapping, decoding and channel estimation using 
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an estimation-maximisation (EM) algorithm and a turbo-demodulation procedure has been proposed. However, the 
proposed algorithm has been derived in a single carrier context and is therefore unsuitable for OFDM implementation. 
[0014] A turbo channel estimation method using EM algorithm in an OFDM context has been proposed in E. Jaffrot, 
M. Siala, "Turbo channel estimation for OFDM systems on highly time and frequency selective channels", Proc. 
5 ICASSP2000. However, the channel estimation method assumes, as a pre-requisite, a phase-shift keyed (PSK) mod- 
ulation scheme. As such, the method is limited in its applications and works only for binary PSK (BPSK) and quadrature 
PSK (QPSK) mapping. Thus, it cannot be applied to high bit rate applications such as HIPERLAN/2 (which uses a 
quadrature amplitude modulation (QAM) scheme). Moreover, it uses pilot carriers, which reduces the application bit 
rate. 

w [0015] Thus, there exists a need in the field of the present invention to provide an improved channel estimation 
method in an OFDM communication system, wherein at least some of the aforementioned disadvantages may be 
alleviated. 

Statement of Invention 

15 

[0016] The present invention provides a method of channel estimation and a receiver for performing channel esti- 
mation as described in the accompanying claims. 

Brief Description of the Drawings 

20 

[0017] 

FIG. 1 illustrates a graph of mean square error (MSE) for a BRAN A channel with a mobile terminal (MT) moving 
at 3m/s; and 

25 

FIG. 2 shows a simplified processing block diagram of a bit-interleaved and convolutionally coded OFDM commu- 
nication system. 

[0018] Exemplary embodiments of the present invention will now be described, with reference to the accompanying 
30 drawings, in which: 

FIG. 3 illustrates a wireless communication unit adapted in accordance with a preferred embodiment of the present 
invention; 

35 FIG. 4 illustrates a channel estimation model applicable to the preferred embodiment of the present invention. 

FIG. 5 illustrates a block diagram of a decoding unit adapted in accordance with a preferred embodiment of the 
present invention; 

40 FIG. 6 illustrates two examples of mapping that can be used in the preferred embodiment of the present invention ; 

FIG. 7 illustrates a model of an OFDM system used in the preferred embodiment of the present invention; 

FIG. 8 illustrates a block diagram of a decoding unit's forward-backward algorithm in accordance with a preferred 
45 embodiment of the present invention; and 

FIG. 9 shows bit-error-rate simulation results of the channel estimator according to the preferred embodiment of 
the invention. 

50 Description of Preferred Embodiments 

[0019] In summary, the preferred embodiments of the invention described below perform blind channel estimation 
for each OFDM symbol, which allows the receiver to track the channel variations, and consequently to increase the 
time between training sequences. In particular, the channel estimate of the frequency channel coefficients for each 
55 OFDM symbol is performed using an Expectation-Maximization (EM) algorithm that includes estimated channel coef- 
ficients using at least one other, and preferably all the other sub-carriers. 

[0020] Referring now to FIG. 3, there is shown, by way of example, a block diagram of a wireless communication 
unit, adapted to support the inventive concepts of the preferred embodiments of the present invention. The wireless 
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communication unit 300 contains an antenna 302 preferably coupled to an antenna switch 304 that provides isolation 
between a receiver and a transmitter chain within the wireless communication unit 300. 

[0021] The receiver chain shown includes receiver front-end circuitry 306 (effectively providing reception, filtering 
and intermediate or base-band frequency conversion). The front-end circuit 306 is serially coupled to a signal process- 
ing function 308. An output from the signal processing function is provided to a suitable output device 310. 
[0022] The signal processing function 308 performs all signal processing functions for the wireless communication 
unit 300, including demodulation, de-mapping, bit deinterleaving, channel estimation and decoding. In accordance with 
the preferred embodiments of the present invention, the signal processing unit 308 has been adapted to perform EM- 
based channel estimation, using probabilities on coded bits. These probabilities may be provided by any channel 
decoding scheme yielding soft outputs, e.g. a so-called BCJR algorithm - see L. Bahl, J. Cocke, F. Jelinek and J. Raviv, 
"Optimal Decoding of Linear Codes for Minimizing Symbol Error Rate", IEEE transaction on Information Theory, March 
1 974. The preferred embodiment uses a turbo-demodulation process, such that it is able to jointly perform de-mapping 
and decoding of the received signal. The signal processing function 308 is further adapted inasmuch as it employs a 
new channel estimation method, as is further described with regard to FIG. 5. 

[0023] The receiver chain shown also includes received signal strength indicator (RSSI) circuitry 31 2 (shown coupled 
to the receiver front-end 306, although the RSSI circuitry 312 could be located elsewhere within the receiver chain). 
The RSSI circuitry is coupled to a controller 314 for maintaining overall subscriber unit control. The controller 314 is 
also coupled to the receiver front-end circuitry 306 and the signal processing function 308 (generally realised by a 
digital signal processor (DSP)). The controller 314 may therefore receive bit error rate (BER) or frame error rate (FER) 
data from recovered information. The controller 31 4 is coupled to the memory device 31 6 for storing operating regimes, 
such as decoding/encoding functions and the like. A timer 318 is typically coupled to the controller 314 to control the 
timing of operations (transmission or reception of time-dependent signals) within the wireless communication unit 300. 
[0024] As regards the transmit chain, this essentially includes an input device coupled in series via transmit signal 
processor 328 to a transmitter/ modulation circuit 322. Thereafter, any transmit signal is passed through a power am- 
plifier 324 to be radiated from the antenna 302. The transmitter/ modulation circuitry 322 and the power amplifier 324 
are operationally responsive to the controller, with an output from the power amplifier coupled to the duplex filter or 
circulator 304. The transmitter/ modulation circuitry 322 and receiver front-end circuitry 306 comprise frequency up- 
conversion and frequency down-conversion functions (not shown). 

[0025] Of course, the various components within the wireless communication unit 300 can be arranged in any suitable 
functional topology that is able to utilise the inventive concepts of the present invention. Furthermore, the various 
components within the wireless communication unit 300 can be realised in discrete or integrated component form, with 
an ultimate structure therefore being merely an arbitrary selection. The preferred decoding and channel estimation 
function can be implemented in software, firmware or hardware, with the function being implemented in a software 
processor (or indeed a digital signal processor (DSP)) that performs the decoding function, merely a preferred option. 
[0026] The transmitter signal processor 328, which is shown as being distinct from the receiver signal processor 308 
for clarity purposes only, includes a convolutional encoder or a turbo-coding function. 

[0027] In accordance with a preferred embodiment of the present invention, an improved channel estimation algo- 
rithm for the OFDM system is included. The channel estimation uses an Expectation-Maximisation (EM) algorithm. 
This is a parametric estimation algorithm, which attempts to estimate the parameters '9' 404 of a system 400 from the 
observations of its output y 406 blindly without knowing its input x 402, as shown in the model in FIG. 4. 
[0028] The EM algorithm is an iterative procedure, which maximizes the likelihood function in two steps: 

(i) A first Expectation step (E-step): which comprises a computation of: 

EJP(x\ y,e (p) )] [1] 

and 

(ii) A second Maximization step (M-step): 



0<#h-i) = argmaxP()/l0,jc w ) [2] 

0 

In a practical implementation, the maximisation (M-)step is performed by maximizing the auxiliary function: 
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Q(9,9 (p) ) = E x [\ogP(x, y|6)|y,9 (p) ] [3] 

[0029] In accordance with the preferred embodiment of the present invention, the auxiliary function of the channel 
estimation algorithm has been adapted to take into account OFDM specificities. Indeed, the auxiliary function has been 
rewritten as: 



Q(H m , H™ ) = E Xm [log P(x m ,y m ,H im '\H m )\y m ,H^} [ 4 ] 

where: stands for the vector /-/with aOon the /77 th component. This allows, when estimating the channel coefficient 
of the m-th sub-carrier, to use and take advantage of the prior estimations of the channel coefficients of the other 

15 frequency domain sub-carriers. 

[0030] In addition, the Expectation step has been configured to be an estimation of the probability on decoded bits, 
which are needed for the M-step. These probabilities are provided by any decoding scheme yielding soft outputs. 
Furthermore, the preferred embodiment of the present invention, proposes to use a turbo-demodulation procedure, in 
order to optimise jointly the de-mapping and the decoding steps. In this manner, the decoding algorithm leads to a 

20 much better performance than a non-iterative decoding scheme. 

[0031] However, it is within the contemplation of the invention that the inventive concepts herein described are ap- 
plicable to any convolutional encoder, i.e. it is applicable for any rate and any puncturing scheme. Furthermore, the 
decoding scheme of the present invention can be applied to any encoded data that can be decoded by a soft output 
decoder, i.e. any convolutional or turbo decoder. 

25 [0032] Referring now to FIG. 5, a block diagram model 500 of the preferred channel estimation procedure is illus- 
trated. The received data stream 505 is demodulated in OFDM demodulation function 510. The output from the de- 
modulation function 510 is input to a channel estimation function, which incorporates an estimation function 520 and 
a maximization function 550. 

[0033] The OFDM demodulated signal is input to a de-mapping function 525, which outputs de-mapped signals to 
30 a deinterleaving function 530, and thereafter to a soft output decoder function 535. In this manner, the channel esti- 
mation function is operational over both the de-mapping and decoding functions. 

[0034] The de-mapping function in the preferred embodiment, i.e. the function that associates a subset of bits to a 
modulation constellation location, is able to operate with any labelling map. For example, as shown in FIG. 6, the de- 
mapping operation may be performed on a subset of a 16-QAM signal, which is composed of four bits. Each subset 

35 has to be associated to one of the sixteen positions of the constellation. There are a number of different ways to define 
the association function, which affect the performance of the decoding system. For example, constellation 61 0 in FIG. 
6 illustrates a Gray labelling structure, which is deemed optimum when a Viterbi decoder is used. However, when turbo- 
demodulation is used, an alternative mapping format such as the set partitioning format 620 may provide better results. 
[0035] The output from the decoding function 535, is fed back to the de-mapping function 525 via an interleaving 

40 function 540. The output from the improved estimation function 520 is input to a maximization function 550. The output 
from the maximization function 550 is input to the de-mapping function 525. 

[0036] It is noteworthy that this arrangement also yields decoded bits, in a bit-stream 535. The bit decisions are made 
from decision function 560. Thus, the improved channel estimation arrangement is a joint channel estimation and 
decoding scheme. 

45 [0037] In accordance with the preferred embodiment of the present invention, the channel estimation operation has 
been adapted to take into account the OFDM specificities, using the coded bits on neighbouring frequency channels 
[0038] Referring now to FIG. 7, a block diagram of a parallel equivalent representation of a cyclic-prefix OFDM 
modulator-demodulator 700 is illustrated. In effect, an OFDM transmitter wishes to emit a set of bits 
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<m<M 



702 over a channel 
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using, say, a 2 N -QAM scheme 704. 
[0039] First every subset of Nbits 



F« Ji< 



}\<n<N 

702 is modulated into a symbol x m . In the OFDM context, the 



\. X m \< 



symbols 706 are transformed from a frequency domain representation into a time domain using an inverse fast fourier 
transform (IFFT). The time-domain signals are then serialized and convolved by the channel filter. The signals are 
effectively corrupted by noise 



\<m<M 



732 before they are finally returned into the frequency domain using a FFT (not shown) to produce an output signal 



{y m \ 



734. 

[0040] It is noteworthy that in a standard OFDM transmitted signal, the signal is confined into a spectrum mask by 
precluding some carriers from being used, especially the ones at the top and the bottom of the FFT. In this regard, 
such carrier values are set to zero. 

[0041] An important characteristic of OFDM communication systems, utilised in the preferred embodiment of the 
present invention, is that OFDM systems are dimensioned so that the length of the channel impulse response ('CIR') 
is smaller than the cyclic prefix size. Accordingly, it is possible to use the classical OFDM system model representation 
of FIG. 7. 

[0042] Let us consider the case where h is a vector of M independent variables, with its first L components being L 
independent Rayleigh fading values with variance of, 1 = 1, L Let us further assume thatthe M-L other components 
are independent Rayleigh fadings with a small variance e 2 . 

[0043] If the variance e is small enough, h has been determined as providing a reasonably good model of the time 
domain channel. This leaves the frequency domain channel as being represented by: 



H = Fh , with F = 



4m 



M I 



[5] 



[0044] In accordance with the preferred embodiment of the present invention, let us define a new auxiliary function 
that provides an improved estimate of a channel coefficient H m . It is noteworthy that the improved channel estimate 
considers the channel coefficients on all other carriers in order to refine/track the primary channel coefficient. 
[0045] The inventors have appreciated that knowledge of the channel coefficients from other carriers is useful, as 
the decoder needs an estimate of the frequency channel coefficients 



fa.}, 



\<m<M 



to eradicate the channel effects from the received coded bits. Indeed, in the preferred embodiment of the present 
invention, the frequency channel coefficients from the Fourier transform of the time channel coefficients 
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i h lhl<L 

are utilised in this regard. 
[0046] Thus: 

H = Fh [6] 

Where: 

H is a column vector of 



l/-*mJl<m<M ' 

h is a column vector of 

{ h llsl<L 

filled with zeros for the M-L last coefficients, and 
Fisa 'MxM' matrix filled with Fourier coefficients. 

[0047] At the opposite end of the time domain, the values of elements of H, are unknown. Indeed, due to the intro- 
duction of white noise in the model, it can be readily assumed that all of the elements are corrupted. 
[0048] The inventors of the present invention have appreciated that a standard channel estimation method, such as 

A y 

H m =-^-, where x m is a known symbol, leads to a poor result. 

x m 

[0049] In contrast, the proposed method of using the knowledge of the channel coefficients from other carriers helps 
to refine the values of channel estimate. In particular, the definition of the statistical model h and the cost function are 
key distinguishing features of the inventive concepts of the present invention. 
[0050] Thus, in the auxiliary function, we add/7 fm> in the following manner: 



Q{H m , tf « ) = E Xm [log nx m ,y m M (m) \HJ\y m ,H { f] [ 7 ] 

[0051] Where: 

«w> stands for the vector Hwith a Con the m th component. 

[0052] The maximization of this function leads, after some computational effort, to the following channel estimate 
update formula: 



[8] 



where: 
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7 2 =Ws H [ m , m) , A SH =SF*T>iag([(T?.--<rle 2 -e 2 ])FS T [9] 



and 



jr<r? + (M-L)eX . [10] 



S is the information carriers selection matrix. When dealing with a full sub-carrier OFDM symbol, 2 can be rewritten into: 

15 

__ 2 1 ^ 2 M - L 2 

M i~f M 

20 

[0053] It is noteworthy that the general expression of y 2 of equation [9] takes into account the fact that some sub- 
carriers are not used, which is always the case in practice. 

[0054] In practice, three parameters are needed to implement the aforementioned formula, namely: P(y m \x m , hfi^), 
P(x m ) and no,- . 

25 [0055] If we assume that the channel additive noise is a Gaussian white noise of variance a 2 P(y m \x m ,hf p) )may be 
viewed as: 

P(yJ*„X P) )~e ° 2 [12] 
no 

[0056] Where: 

35 P( x m) ' s th e product of the probabilities of each coded bit involved in the symbol x m . The turbo-demodulation 

process provides an estimation of these probabilities. Thus, P(x m ) is the result of the E-step; and 

is unknown, but each of its component may be taken as the current estimate of the channel coefficient, i. 
e. nr = H (p+7) for i<m and fi« = H (p) for i>m. 

40 

[0057] However, the inventors of the present invention have recognised that the above methodology would lead to 
the last set of coefficients being estimated more accurately than the first set of coefficients. Hence, to remove this 
drawback and improve the estimation, the inventors propose implementing a 'forward-backward' procedure to further 
improve the channel refinement. 

45 [0058] The preferred implementation of the forward/backward process is illustrated in FIG. 8. The channel refinement 
is performed on a carrier-by-carrier basis. Received data is input to the Turbo demodulation estimation process 804, 
as illustrated in FIG. 6. Furthermore, the received data is also input to the maximization process 806, which comprises 
an algorithm having two steps, i.e. one 'forward' step 808 and one 'backward' step 810. The operation of the forward 
step 808 and the backward step 810 is performed on the received data and the channel estimation values output 812 

50 from the estimation algorithm 804. The operation of the forward step 808 and the backward step 810 is performed 
making bit decisions 814 using the information bit probabilities. The output of the process is a sequence of received 
bits 816. 

[0059] In accordance with the preferred embodiment of the present invention, the forward step is implemented as 
follows: 

55 The H m coefficients are estimated from m=1 to m=M. 

[0060] In this step, when computing the channel coefficient on sub-carrier m, the channel coefficient estimates on 
the other sub-carriers 1 , such that m<1<M, come from the previous iteration. The channel coefficient estimates on the 
other sub-carriers 1 such that 1 ('one') <1<m come from the actual iteration. 



10 



v 2 =- 



M 



8 



EP 1 416 688 A1 



[0061] The backward step is implemented as follows: 

The H m coefficients are estimated from m=M\.o m=1, using the estimates of the forward step. 
[0062] In this step, when computing the channel coefficient on sub-carrier m, the channel coefficient estimates on 
the other sub-carriers 1 such that m<1<M come from the forward step, and the channel coefficient estimates on the 
5 other sub-carriers 1 such that 1 ('one') <1<m come from the present backward step. 

[0063] In this manner, approximately the same level of accuracy of each channel coefficient can be achieved, thereby 
improving the channel estimate performance. 

[0064] To fully evaluate the performance of this method, a set of Monte-Carlo simulations was performed over two 
hundred channels to estimate the bit error rate (BER) and the packet error rate (PER). The simulation context was 
10 performed using one hundred OFDM symbols of sixty-four sub-carriers preceded by two pilot symbols. This is useful 
in obtaining a first channel estimate. The simulation was performed over sixteen independent Rayleigh fading channels. 
A half rate convolutional encoder (035, 023) was used. Symbols were bit-interleaved with a random pattern and mod- 
ulated by a 16-QAM constellation. 

[0065] The results of the simulation are illustrated in FIG. 9 and FIG. 10. Referring now to FIG. 9, a comparison of 
15 the average BER versus signal to noise ratio (SNR) is illustrated for an OFDM system. FIG. 1 0 illustrates a comparison 
of the average PER versus SNR for the OFDM system. 

[0066] In both FIG. 9 and FIG. 10, the OFDM system does not refine channel estimates and performs only two 
iterations of the turbo-demodulator per global iteration. Three curves are shown: 

20 (j) A first curve indicates the performance without employing the inventive concepts herein before described, iden- 

tified as 'NONE 2 TD'; 

(ii) A second curve illustrates the performance when implementing the proposed invention, identified as '1 
EM-OFDM-SS 2 TD for e 2 =3e" 3 ; and 

(iii) A third curve illustrates the performance when implementing the proposed invention, identified as '1 
25 EM-OFDM-SS2 2 TD' for e 2 =1 e" 1 6 . 

[0067] It is noteworthy that, by implementing the inventive concepts herein described, an improvement of around 
2dB in BER/PER can be achieved. Moreover, provided that the variance e parameter is small enough the setting of 
the parameter is not critical. Hence, the user does not have to perform precise tuning on the variance e parameter. 
30 [0068] It is within the contemplation of the invention that the aforementioned inventive concepts can be applied to 
any element in the communication system that performs decoding and channel estimation, for example a BTS com- 
munication unit and/or a MS. 

[0069] In the preferred embodiment of the present invention, the maximization step utilises a forward-backward ap- 
proach. However, it is within the contemplation of the invention that alternative maximization techniques may be used 
35 in the decoder such that the decoder is still able to benefit from the inventive concepts described herein. 

[0070] It will be understood that a method of channel estimation performed by a communication unit operating in an 
OFDM communication system as described above, tends to provide at least some of the following advantages: 

(i) A mechanism to perform jointly de-mapping, decoding and channel estimation is described, whereby their per- 
40 formancecan be globally optimised. In particular, the modification of the EM cost function results in the modification 

of the channel update formula, which provides an improved decoding performance. 

(ii) There is no estimation error in a noiseless case. In practice, this means that the proposed estimator performs 
better than the strongest-path method based estimator for high SNR conditions. 

45 

(iii) Estimation noise is partially removed from the selected channel coefficients. 

(iv) There is minimal additional complexity when compared to prior art arrangements. 

50 (v) The inclusion of OFDM specificities in a novel and inventive approach enables bit-interleaved coded modulation 

improvements proposed by J. Boutros et al. to be translated to an OFDM scenario. In particular, the inventive 
concepts can be applied to any kind of mapping, and does not use pilot carriers. Nevertheless, the inventive 
concepts hereinbefore described may be used in the presence of pilot carriers, if desired. 

55 [0071] Whilst specific, and preferred, implementations of the present invention are described above, it is clear that 
one skilled in the art could readily apply further variations and modifications of such inventive concepts. 
[0072] Thus a method of channel estimation performed by a communication unit operating in an OFDM communi- 
cation system has been provided wherein at least some of the aforementioned disadvantages with prior art arrange- 
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merits have been alleviated. 



Claims 

1. A method of channel estimation in a wireless orthogonal frequency division multiplexed (OFDM) communication 
system (700), comprises the steps of: 

receiving a signal in the time domain; 

applying a Fourier transform to said received signal to obtain a frequency domain signal including a plurality 
of sub-carriers; 

estimating probabilities of coded bits for at least said plurality of frequency domain sub-carriers; and 
performing channel coefficient estimation for at least said plurality of frequency domain sub-carriers 

characterised in that 

said step of performing channel coefficient estimation for each of said plurality of frequency domain sub- 
carriers uses channel coefficient estimates for at least one other of said plurality of frequency domain sub-carriers. 

2. A method of channel estimation according to Claim 1 , wherein said step of performing channel coefficient estimation 
for substantially each of said plurality of frequency domain sub-carriers uses channel coefficient estimation benefits 
from said channel coefficient estimates for substantially all the other frequency domain sub-carriers of said plurality. 

3. A method of channel estimation according to Claim 2, wherein said plurality of frequency domain sub-carriers 
comprises substantially all the sub-carriers of said frequency domain signal. 

4. A method of channel estimation according to any preceding Claim and comprising repeating said steps of estimat- 
ing probabilities and performing channel coefficient estimation so as to improve iteratively an accuracy of said 
channel coefficient estimates. 

5. A method of channel estimation according to Claim 4, wherein the kth channel coefficient estimation is substantially 
in accordance with the following equation: 



where is the (p+7)th estimate and the pth estimate of the channel coefficients, y k is the received data 

corresponding to the transmitted data x k , a 2 is the channel noise variance,^ is the channel coefficient vector H 
with a 0 on the /c th component and A -1 , u 2 and y 2 have the meanings indicated hereinabove. 

6. A method of channel estimation according to any of Claim 4 or 5, wherein the step of performing channel coefficient 
estimates comprises replacing previously estimated channel coefficients of said plurality of frequency domain sub- 
carriers with respective current channel coefficient estimates. 

7. A method of channel estimation according to any of Claims 4 to 6, wherein repeating said step of performing 
channel coefficient estimation comprises applying a cost function on an Expectation-Maximization algorithm on 
said plurality of frequency domain sub-carriers to improve said channel coefficient estimates. 

8. A method of channel estimation according to Claim 7, wherein said step of performing a channel coefficient esti- 
mation includes calculating an auxiliary function, the method further characterised by the step of: 

performing a Maximisation process on said auxiliary function in substantially the following manner: 
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Q(H m ,HL P) ) = Ejlog?(x m ,y m ,H™ \H m ) I y m ,H<»]. 

A method of channel estimation according to any of Claims 4 to 8, wherein said step of performing a channel 
coefficient estimation includes applying a forward -backward algorithm on said received signal to said plurality of 
channel coefficient estimates in which estimates are made in a first order of said plurality of frequency domain 
sub-carriers and subsequently estimates are made in a reversed order of said plurality of frequency domain sub- 
carriers so as substantially to equalise an estimation accuracy across said plurality of frequency domain sub- 
carriers. 

0. A method of channel estimation according to any preceding Claim, wherein said step of channel estimation is 
performed on received data that is decoded by a soft output decoder, said decoder providing estimations of said 
probabilities. 

1. An orthogonal frequency division multiplexed (OFDM) receiver for a method of channel estimation as claimed in 
any preceding Claim, and comprising: 

demodulation means for applying said Fourier transform to said received signal to obtain said frequency do- 
main signal including a plurality of sub-carriers; 

decoding means for decoding the received signal and estimating said probabilities of coded bits for at least 
said plurality of frequency domain sub-carriers; and 

channel estimation means for performing channel coefficient estimation for each of said plurality of frequency 
domain sub-carriers using channel coefficient estimates for at least one other of said plurality of frequency 
domain sub-carriers. 
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